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Abstract 

We discuss the first step in the moduli stabilization program a la KKLT for a 
general class of resolved toroidal type IIB orientifolds. In particular, we discuss 
their geometry, the topology of the divisors relevant for the D3-brane instantons 
which can contribute to the superpotential, and some non-trivial aspects of the 
orientifold action. 


1 Introduction 

Most string theory compactifications come with a large number of moduli. Corresponding 
to massless helds in the low energy effective action, they give rise to a hfth force of about 
gravitational strength and are consequently in conflict with experiment. Therefore, a lot 
of effort is being dedicated to the subject of moduli stabilization. 

Compactihcations which should ultimately yield the MSSM are required to have N = 

I supersymmetry in four dimensions. One way to achieve this is to start with type 

II string theory compactihed on a compact CalabmYau threefold X yielding N = 2 
supersymmetry. Having non-abelian gauge groups requires the addition of D-branes 
wrapping supersymmetric cycles in X. Their presence induces an RR flux which has to 
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be cancelled by adding appropriate orientifold planes. Withont D-branes and orientifold 
planes, the relevant modnli are the dilaton, and the deformations of the complex structure 
and the Kahler form of X while the moduli of the RR helds decouple. After adding the 
orientifold planes, the geometric and the RR moduli get mixed non-trivially. 

In the KKLT-proposal ^ , starting from type IIB theory, the dilaton and the complex 
structure moduli are stabilized via background three-form fluxes, while the Kahler moduli 
are hxed by non-perturbative effects. In this way, a stable, supersymmetric AdS vacuum 
can be obtained. In a second step, this vacuum is supposed to be lifted up to a metastable 
deSitter vacuum by introducing anti-D3-branes, which give rise to a positive contribution 
to the potential, producing the uplift to a small value of the cosmological constant. 

We will be concerned with the hrst step only, i.e. we study the issue of the stabilization 
of all moduli. We will present a program for testing this proposal in a number of examples 
and point out issues which deserve to be understood more deeply in order to fully realize 
it. This is a continuation of [2] and a preview of some of the results of jS]. 

The stabilization of the complex structure moduli and the dilaton via three-form fluxes 
has been extensively discussed in the literature |3] (for a review, see jS] and references 
therein), while the general conditions for the generation of a non-perturbative super¬ 
potential depending on the Kahler moduli have recently received a lot of attention, 
see iniiiiiHiiniiiniiiiiiiaiiniiii!. On the other hand, some concrete realizations of 
step one of KKLT, i.e. models with all moduli hxed have been put forward nn, 

ini. 

The closed string moduli space consists of Kahler moduli, which parametrize 

the size of X, h‘^’^{X) complex structure moduli, parametrizing the shape of X, plus the 
complex dilaton S = i Cq + e~'^, which consists of the RR 0-form Cq and the dilaton <p. 
The three-form hux that is turned on is of the form + iS H^, with = dC 2 is 

the held strength of the RR 2-form C 2 and = di ?2 is the held strength of the B-held. 
Such a 3-form hux gives rise to a superpotential ^Hj, IIH] of the form 

hTflux= [ GsAfi, (1.1) 

Jx 

where is the Calabi-Yau (3,0)-form. This superpotential depends on the dilaton 
through the hux, and on the complex structure moduli through but since there is 
no dependence on the Kahler moduli, one must look for other means to stabilize them. 
Here, non-perturbative ehects yield such a means. There are two possible origins of such 
a non-perturbative superpotential: Euclidean D3-brane instantons wrapping 4-cycles 
D G X and gaugino condensation in the world-volume of D7-branes wrapping again 
4-cycles D G X. Both give rise to terms in the superpotential of the form 

( 1 . 2 ) 
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where Vi is the volume of the divisor Di wrapped by the brane and therefore depends on 
the Kahler moduli. The prefactor Qi comes from a one-loop determinant and in general 
depends on the complex structure moduli. Unfortunately, it is very hard to calculate 
explicitly, however recently a nice method has been presented in HU Otherwise, the 
best one can do in general is deciding whether it is zero or not. In the absence of fluxes 
Witten gave a criterion for this in the framework of F-theory The Calabi-Yau 

threefold X is lifted to a Calabi-Yau fourfold Y and the divisor D (Z X is lifted to 
divisor D C Y. Then the condition is that the holomorphic Euler characteristic of D has 
to equal one: 

\'(Ob)=Y'»°'‘(F = 1. (1.3) 

i=0 

Under favorable circumstances, this condition can be reduced to a condition on the 
holomorphic Euler characteristic of the divisor D. The Hodge numbers and 

are also relevant to decide whether there are contributions from D7-branes |2U. 


2 Moduli stabilization for resolved toroidal orbifold 
models 

We want to know whether the KKLT scenario can actually be realized, and in a hrst 
step, one is trying to understand the problem of moduli stabilization. One way to gain a 
better understanding is to construct specihc examples which £t all requirements. While 
the complex structure moduli and the dilaton can be generically hxed by fluxes, the 
hrst impression was that divisors which generate contributions to the non-perturbative 
superpotential are not generic [20], [22] • 

In the search for models with all moduli fixed by the effects described above, we 
follow the example of and consider orientifolds of resolved toroidal orbifolds. The 
motivation to look at this specific class of models is that toroidal orbifolds have a simple 
structure which allows for many quantities to be computed explicitly. Furthermore, they 
are also phenomenologically attractive because of their standard model-like features such 
as non-abelian gauge groups, chiral fermions and family repetition. 

Orbifolds come with a number of untwisted Kahler moduli, which are inherited directly 
from the underlying torus, plus a number of twisted Kahler moduli which parametrize 
the sizes of the exceptional divisors which result from the resolution of the orbifold 
singularities. If one takes the aim of hxing all moduli seriously, one has to take care 
of the twisted Kahler moduli as well. Most likely, they will not be stabilized at zero 
size, so one has to look at the smooth Calabi-Yau manifold resulting from all orbifold 
singularities being resolved. The resolution process is most conveniently described in the 
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language of toric geometry. Divisors originating from blow-ups are particularly promising 
candidates for yielding a non-perturbative superpotential since they are in general rigid, 
i.e. have no moduli. 

Ultimately, we would like to check all toroidal orbifolds which allow tadpole cancellation 
in type IIB for their suitability for the KKLT proposal. The example of T®/Z 2 x Z 2 which 
was studied in great detail in [ini turned out to allow stabilization of all moduli. In the 
class of toroidal orbifolds, this model holds a special place because of its simplicity (the 
orbifold group contains only three non-trivial elements, which give rise to only one sort 
of local patch) and symmetry (the conhguration of hxed sets as well as of the 0-planes is 
the same for every coordinate direction) and its F-theory lift is known. This is radically 
different for most other toroidal orbifolds as we will explain in the following sections. For 
this reason we would like to ascertain how common the stabilization of all moduli in the 
class of toroidal orbifolds is. Furthermore, we believe that the construction of blowing 
up the various different patches and gluing them together to form a smooth Calabi-Yau 
manifold is worthwhile to be understood in its own right. 

In the following, we will sketch the recipe of and go on to give some of the steps 
explicitly for our example. 

• We begin with some toroidal orbifold. The regions close to the orbifold singularities 
can be described with toric geometry. In these local patches, we can resolve the 
singularities via a blow-up. 

• In the next step, those local patches have to be glued together to form a smooth, 
compact Calabi-Yau threefold X. 

• The topologies of the divisors of X have to be analyzed. This will be important in 
order to decide which of them can generate a non-perturbative contribution to the 
superpotential. 

• A consistent orientifold projection has to be found and performed. This yields 
the 0-planes and changes the geometry. The tadpoles can be cancelled by adding 
D-branes. 

• The Kahler potential has to be determined. For the Kahler part, we have to 
compute the intersection ring of X. This yields a non-trivial consistency check on 
the the previous steps. 

• The 3-form flux is turned on to £x the complex structure moduli and the dilaton. 
The effect of the background flux on the fermionic zero modes on the D-branes 
and therefore the generation of the non-perturbative terms in the superpotential 
has to be taken into account. 
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Group element 

Order 

Fixed set 

Conj. classes 

9 

6 

12 hxed points 

12 

02 

3 

9 hxed lines 

6 

03 

2 

16 hxed lines 

8 


Table 1: Fixed point set for Zq_jj on SU{2Y x SU{3) x G 2 


• In the last step, the scalar potential obtained from the superpotential and the 
Kahler potential has to be minimized. 


2.1 Toroidal Orbifolds 


Probably the simplest way to construct a Calabi-Yau compactihcation is to take a 
toroidal orbifold and resolve its singularities. Such orbifolds have been analyzed in detail 
some time ago. This analysis resulted in a partial classification 1231, m which is our 
starting point. These orbifolds are of the form T^/G where G is either or x 

where G acts by multiplying the complex coordinates by phases. The G action has fixed 
sets at which the orbifold is singular. A point is hxed under 6"^ ^ G, n = 0,..., |G| — 1, 
if it fulfills 


«“/<"> = /<"> +a, oeA, 


( 2 . 1 ) 


where a is a vector of the torus lattice A. These singularities have to be resolved in order 
to obtain a smooth Calabi-Yau threefold and we review this in the next subsection. 

As a specific example, we will discuss here T^/Zg-// compactified on the root lattice 
of SU{2Y X SU{3) X G 2 - In complex coordinates, the twist acts as 

9: ^ {e e = (2.2) 


The hxed point sets are listed in Table 12.11 and a schematic picture is displayed in 
Figure n Each complex coordinate is shown as a coordinate axis and the opposite faces 
of the resulting cube are identihed. The covering space is shown, not the quotient. 

Here we see a hrst important difference from the example in ng. The conhguration 
of the hxed sets is diherent in the diherent coordinate directions and there are several 
diherent types of singularities. This is the typical situation for all the orbifold models 
we consider. Furthermore, it is not yet known how the F-theory lift can be explicitly 
constructed. From the hxed point conhgurations one hnds the Hodge numbers to be 
= 3 + 32 = 35 and = 1 + 10 = 11. From these one can see that this manifold 
possibly allows for a Borcea-Voisin realization Bn.ra, from which the F-theory fourfold 
could be constructed. This is, however, not the case in general. If no F-theory lift is 
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Figure 1: Schematic picture of the hxed set conhguration of Ze_ II on SUi2yxSU{3)xG2 


known, Witten’s criterion (El cannot be applied. A possible way out is the method 
presented in [T^ . 

2.2 Resolution of Singularities 

A systematic way to resolve abelian orbifold singularities is toric geometry [2^]. In the 
following paragraphs, we summarize some of the basic facts about toric geometry. An 
n-dimensional toric variety takes the form 

Xs = (C"\Fs)/(e)^ (2.3) 

where n = d — r, and (C*)'’ acts by coordinatewise multiplication. The set F-£ is a 
subset that remains hxed under a continuous subgroup of (C*)'' and must be subtracted 
for the variety to be well dehned. The action of (C*)^ is encoded in a lattice N which 
is isomorphic to Z'’* and by its fan S, a collection of strongly convex rational cones in 
N IR with the property that each face of a cone in S is also a cone in S and the 
intersection of two cones in S is a face of each. The fc-dimensional cones in S are in 
one-to-one correspondence with the codimension fc-submanifolds of X^. In particular, 
the one-dimensional cones correspond to the divisors in Xs- The fan S can be encoded 
by the generators of its edges or one-dimensional cones, i.e. by vectors Vi G X. To each 
Vi we associate a homogeneous coordinate Ui of The (C*)'’ action is encoded on the 
Vi in r linear relations 

d 

= a = l,...,r, (2.4) 

i=l 
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We are only interested in Calabi-Yan orbifolds C^/G of dimensions m = 2, 3, so we 
reqnire to have trivial canonical class. This translates to demanding that all bnt one 
of the Vi lie in the same affine hyperplane one nnit away from the origin wq- This means 
that the last component of all the Vi (except no) eqnals one. This allows ns to draw toric 
diagrams in two (one) dimensions instead of m = 3 (m = 2). The fan E associated to the 
former is obtained as follows: We have a single three-dimensional cone in S, generated 
by fi, ^ 2 , ns. A generator 6* of G of order n acts on the coordinates of by 

9: (2.5) 

Then the local coordinates of Xy, are (^^3Yv3)k_ coordi¬ 

nates in N of the generators n* of the fan, we reqnire the to be invariant nnder the 
action of 9. This resnlts in finding two linearly independent solntions of the eqnation 

{vi)k + ni {v 2 )k + n 2 {vYk = 0 modn. (2.6) 

The divisors corresponding to Vi will be denoted by Di. 

Xy is smooth if all the top-dimensional cones in S have volnme one. Here, there is only 
one snch cone whose volnme is |G|, hence Xy is singnlar. There is a standard procednre 
for resolving singnlarities of toric varieties. It consists of adding all lattice points in N 
which lie in the polyhedron in the affine hyperplane at distance one which is spanned 
by the generators Vi. For the fan E this means that the corresponding one-dimensional 
generators Wi are added to it and that it has to be snbdivided accordingly. We denote 
the refined fan by E. In it is shown that these new generators can be related (in the 
case m = 3) to the generators of G as follows: 

Wi = V2 + gf ns, gt'^ = 1, 0 < g^Y^ < 1. (2.7) 

k=l 

where = {gi \ <72*^ S's*^) ^ represents the corresponding generator 9i. The cor¬ 

responding exceptional divisors are denoted Ei. The snbdivision of the fan E into E 
corresponds to a triangnlation of the toric diagram. In general, there are several trian- 
gnlations, and therefore several possible resolntions. They are all related via birational 
transformations. 

The case m = 2 is even simpler. The singnlarity is called a rational donble 

point of type A„_i and its resolntion is called a Hirzebrnch-Jnng sphere tree consisting 
of n — 1 exceptional divisors intersecting themselves according to the Dynkin diagram of 
A„-i. 

In onr example, we see from Table ^ that we have lines of C^/Z 2 and C^/Zs singn¬ 
larities, as well as singnlar points of the form C^/Zg. For the latter we take the action 
in (12.2p . Then (12.6p leads to the following three generators of the fan E: 

m = (-2, -1,1), V2 = (1, -1,1), ns = (0,1,1). (2.8) 


7 


The points inside the triangle {vi,V 2 .,v^) lead to fonr additional generators of S: 

n;i = (0,0,1), n; 2 = (0,-l,l), ^3 = (-1,0,1), n;4 = (-1,-1, !)• (2.9) 


By (12.711 they correspond to the generators 6 **, i = 1,..., 4, of In this example, 

there are five triangulations. Figure El shows three of the corresponding toric diagrams. 



Figure 2: Toric diagrams of the resolutions of C^/Zg-// 

The remaining two not shown here are obtained by (i) taking the second case and blowing 
the curve Ei ■ down and then blowing the curve E 2 ■ D 2 up (such a process is called a 
flop) and (ii) taking the third case and flopping the curve Ei ■ E 2 to ■ E 4 ^. 

2.3 The global resolution 

Having analyzed the local structure of the singularities and their resolution we can return 
to the resolved global orbifold. The general procedure consists of gluing all the resolved 
local patches Xg. This is achieved by relating their non-compact divisors Di to the 
divisors Ri, which are inherited from the covering torus T®, via the linear relations of the 
toric varieties Xg. 

The inherited divisors and the exceptional divisors together form a basis for the divisor 
classes of the resolved orbifold. In addition, there are further natural divisors which will 
be expressed in terms of linear combinations of these basis elements. These additional 
divisors are simply planes lying at fixed loci of the orbifold action: Dia = {z^ = 
where a runs over the fixed loci in the ith direction. In terms of the local toric patches 
they correspond to the non-compact divisors D^. The local linear relations become “slid¬ 
ing” divisors in the compact geometry, and are related to the inherited divisors as follows. 
Consider the divisors {z^ = c 7 ^ ^fixedol- They “slide” in the sense that they can move 
away from the fixed point, and the way they can move is constrained by the linear re¬ 
lations in the local geometry. We need, however, to pay attention whether we use the 
local coordinates 5* near the fixed point on the orbifold or the local coordinates on 
the cover. Locally, the map is 5* = (^*)”, where n is the order of the group element that 
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fixes the point. The divisor Ri = {R = c"'} on the orbifold lifts to a union of n divisors 
R^ = with e” = 1. Consider the local toric patch before blowing up. 

The hxed point lies at c = ^ and in the limit as c approaches this point we hnd the 

relation between Ri and Di to be Ri ~ nDi. This expresses the fact that at the fixed 
point the polynomial defining Ri on the cover has an nth order zero on Di. In the local 
toric patch Ri ~ 0, hence nDi ~ 0. After blowing up, the Ri and nDi differ by the ex¬ 
ceptional divisors Ek introduced in the process of resolution. The difference is expressed 
precisely by the linear relation (Q corresponding to the resolved toric variety Xg and 
takes the form Ri ~ nDi + From this we see that this relation is independent 

from the chosen resolution. Such a relation holds for every fixed point a which adds 
the label a to the relation: 

Ri ~ nDia + Eka- (2T0) 

k 

The precise form of the sum over the exceptional divisors depends on the singularities 
involved. 

For our example, we obtain the following: From Figure ^ we see that we have 12 local 
patches which each sit at the intersection of a C^/Z 2 and a C^/Zs fixed line. 
They each have one exceptional divisor in the interior and three on the boundary of the 
toric diagram. From the discussion in Section f2.2l we see that the 6 C^/Za fixed lines 
each contribute two exceptional divisors. At the C^/Zg-// fixed point sitting on this line, 
they are identified with the two divisors sitting on the boundary {Di,D 2 ). Then, there 
are 8 lines fixed under the Z 2 -element, contributing each one exceptional divisor. In the 

/Xq_jj fixed point sitting on top of such a fixed line, the divisor is identified with the 
one on the boundary {Di, D 3 ) of the toric diagram of the resolution, as before. We have 
illustrated this in Figure El The vertices corresponding to the exceptional divisors are 
colored to match Figure H This gives a total of 12 ■ 1 -|- 6 ■ 2 -|- 8 ■ 1 = 32 exceptional 
divisors. In addition, we have three inherited divisors Ri = {5* = c®}, z = 1,2,3 which 
yields = 35. 

Again, from Figure [T] we see that the divisors located at fixed sets are three planes 
T* 2,/3 at the Z 3 fixed lines and four planes D^^.y at the Z 2 fixed lines. In the direction 
some of the fixed lines get identified by 6, so that there remain three equivalence classes 
Di a- As representatives we choose 2 ;^ = 0 for a = 1, z^ = 1/2 for a = 2 , and z^ = 1/3 
for a = 3. So for E^^^a, we can have a = 1,2 and for E 2 / 4 ,^f 3 a, we can have a = 1,3. 
Using in dzi, the relations (ITTUl) become 

4 34 3 

= 6 Dl l + 3 -^ 3 , 7 ! + [2 ^2,131 + 4 E^^pi], 

7=1 / 3=1 7=1 13=1 

4 

Ri = 6 Di^2 + 3 .^ 3 , 72 , 

7=1 
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Figure 3: Toric diagrams of the resolution of the hxed point in T^/Zg-//. We 

emphasize that the points in the interior of the boundary are identihed with the Dynkin 
diagrams of Ai and A 2 obtained from resolving the C^/Z 2 and hxed lines which 

intersect in this C^/Ze_// hxed point. 

4 

Rl = 6 ^1^3 + [2 i?2,/33 +-£'4,/33]) (2-11) 


/3=1 
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i?2 = 3 1^2,/3 + -^1,/37 + 

[2 -^2,/3a 

+ -S4,/3a], 

/3 

7=1 

Q 

a=l,3 



O 

i?3 = 2 -Ei,/37 + 

^ E-i 

7 = 1,.. 

.,4. 

/3=1 

a=l,2 




These relations encode the gluing of the local patches and will be crucial in determining 
the intersection ring of the orientifold. 

2.4 Topology of the divisors 

At this point, we can already determine the topology of the divisors in order to decide 
whether they can contribute to the superpotential (after the lift to F-theory). This can 
be done straightforwardly for all the cases. 

The topology of the exceptional divisors can be discussed in the local toric variety Xg. 
For that purpose we recall the notion of the star of a cone a, denoted Star(cr), which 
is the set of all cones r in the fan E containing a. In our situation, the topology of 
an exceptional divisor Ei is then determined in terms of Star((Tu,J. This means that we 
simply remove from the fan S all cones, i.e. points and lines in the toric diagram, which 
do not contain Wi. It can be shown that the compact exceptional divisors (i.e. those in 
the interior of the toric diagram) can either be a P^, a Hirzebruch surface F„, a toric 
blow-up thereof, or a flop of any of these possibilities. For the exceptional divisors on 
the boundary of the diagram the geometry is effectively reduced by one dimension. The 
only compact toric manifold in one dimension is P^, the extra non-compact direction is 
C, and there is again the additional possibility of toric blow-ups and flops. 
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After the gluing process described in Section ESI the non-compact factor is identihed 
with a hxed line 1%^ After the resolution, the topology becomes a blown up at 
some points. Hence, the topology of these divisors is a (possibly non-trivial) product of 
two blown-up P^s. The same argument also applies to the locally non-compact divisors 
Di. In the global orbifold they take the form T^/Z„ x I^Lra- With the singularities 
removed, this is the same as a (possibly non-trivial) product of two P^s with a certain 
number of points removed in each factor. Blowing up the singularities glues in either 
points or further P^s depending on the number of exceptional curves ending in Di in the 
toric diagram. 

The upshot is that all the divisors Di and are rational surfaces which have in 
particular = 0. The inherited divisors Ri naturally come with the position 

modulus c, i.e. have 7 ^ 0 . 

In our example, the only compact exceptional divisor is Ei. In the triangulation b), 
we recognize its star to be that of an Fi. In the other triangulations, Ei is birationally 
equivalent to Fi, since the triangulations are all connected to each other through flops. 
As for the locally non-compact exceptional divisors, in triangulation a) they all have the 
topology of P^ X C, while in the other triangulations, we have hrst to perform blow-ups 
and flops to reduce them to this simple form. 

2.5 Orientifold projection 

In order to obtain a consistent theory, an orientifold projection has to be performed. 
The easiest way would be to start with the corresponding orientifold at the orbifold 
point and then blowing up. However, in this process non-perturbative states from D- 
branes wrapping vanishing cycles may appear, making perturbation theory singular and 
thus preventing the orbifold phase and the smooth phase from being continuously con¬ 
nected (for explicit examples in a slightly different context, see [IHl)- Therefore, 
we have to resolve the singularities hrst and impose the orientifold projection directly 
on the smooth geometry. Since we are interested in the situation with 03- and 07- 
planes, we choose the geometric part of the orientifold action OJg on the torus T® to 
be Jg : (— 2 :^, — 2 :^, —z^). We therefore have to study the hxed sets under 

the combination of Jg and the elements of the orbifold group G. The condition m is 
replaced by 

z = z + a, a G A, i = 0, ...,n — 1. (2-12) 

In contrast to the simple example in uni, these hxed point sets in general do not coin¬ 
cide with the hxed point set of the orbifold group. In particular, there are orbifold hxed 
points which are mapped into each other by the orientifold action. Consequently, the ori¬ 
entifold action is locally trivial. Furthermore, in the split of 
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into an invariant and an anti-invariant part, we now have h}2^ 7 ^ 0. The corresponding 
moduli are not Kahler moduli anymore but have the form ra 

G“ = CJ-SB5. (2.13) 

Currently, there is no superpotential known which could £x these moduli 

In addition, we have to decide whether and how we impose an orientifold action on 
the local coordinates ijk of the exceptional divisors Ek- Presently, no criteria are known 
how to do this except for consistency at the very end of the computations. This leaves 
a lot of freedom. It turns out that indeed a non-trivial action on the i/k is needed. 
In order to emphasize the importance of this point, we mention two places where this 
choice becomes important. When a divisor is hxed by the orientifold action, a factor of 
1/2 appears in front this divisor in ()2.1()|1 resp. ()2.11|1 . Moreover, although the divisors 
of all models in question have similar properties, namely = 0 , it depends 

strongly on the choice of orientifold involution whether such a divisor will contribute to 
the non-perturbative superpotential or not. 


3 Outlook 

Of the seven steps outlined in the beginning of section El we have addressed only the first 
four in this article. The remaining three, such as the calculation of the Kahler potential, 
turning on 3-form flux and the minimization of the scalar potential will be addressed 
in [2], also for the other examples discussed in [2]. 

We can make a preliminary comment on the second but last step of turning on the 
3-form fluxes. Since most of the divisors in our model have vanishing and the 
background flux will have no effect on the zero modes. The reason is that the flux can 
lift zero modes. Since = 0, it can only have an effect on which is 1 

for connected surfaces. But as long as we turn on supersymmetric fluxes, is left 
untouched. 

As a further step a complete formulation of these models in F-theory would be useful 
in order to apply the methods of im to explicitly compute the non-perturbative contri¬ 
butions to superpotential as well as to honestly take into account quantum corrections 
by going away from the constant dilaton limit. The simultaneous simplicity and non¬ 
triviality of these models make them promising candidates for such investigations. 


^We thank Thomas Grimm for pointing this out to us. 
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